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The dark photon, A', and the dark Higgs boson, h !, are hypothetical constituents featured in a 
number of recently proposed Dark Sector Models. Assuming prompt decays of both dark particles, 
we search for their production in the so-called Higgs-strahlung channel, e + e“ —¥ A'h ', with h! —> 
A 1 A 1 . We investigate ten exclusive final-states with A' —> e + e~, , or n + n~ , in the mass ranges 

0.1 GeV/c 2 < my < 3.5 GeV/c 2 and 0.2 GeV/c 2 < my < 10.5 GeV/c 2 . We also investigate three 
inclusive final-states, 2(e + e“)X, 2 (p + and (e + e~)(p, + fj,~)X, where X denotes a dark photon 
candidate detected via missing mass, in the mass ranges 1.1 GeV/c 2 < my < 3.5 GeV/c 2 and 
2.2 GeV/c 2 < my < 10.5 GeV/c 2 . Using the entire 977fb -1 data set collected by Belle, we observe 
no significant signal. We obtain individual and combined 90% credibility level upper limits on the 
branching fraction times the Born cross section, B x (TBom, on the Born cross section, OBorn, and on 
the dark photon coupling to the dark Higgs boson times the kinetic mixing between the Standard 
Model photon and the dark photon, old x e 2 . These limits improve upon and cover wider mass 
ranges than previous experiments. The limits from the final-states 3(7r + 7r _ ) and 2(e + e“)A are the 
first placed by any experiment. For old equal to 1/137, my < 8 GeV/e 2 , and my < 1 GeV/c 2 , we 
exclude values of the mixing parameter, e, above ~ 8 x 10~ 4 . 

PACS numbers: 12.60.-i,14.80.Ec,14.60.-z,14.40.Aq 


Recent results from dedicated dark-matter searches P- 
Sj, muon-spin precession measurements j3], and space- 
based particle observatories Ml may be interpreted as 
deviations from the Standard Model (SM) of particle 
physics. Attempts at devising unified explanations have 
led to Dark Sector Models (DSM) that introduce a new 
hidden or dark U(l) interaction that imbues dark matter 
with a novel charge l8ll28|. A possible mediator of this 
new Abelian force is the dark photon, which has an ex¬ 
pected mass of the order of MeV/c 2 - GeV/c 2 and has a 
very small kinetic mixing with the Standard Model pho¬ 
ton, e, of the order of 10 _5 -10 -2 [15] , The dark U(l) 
symmetry group could be spontaneously broken, by a 
Higgs mechanism, adding a dark Higgs boson h' (or sev¬ 
eral of these) to such models [M] ■ 


Due to the small coupling to SM particles and the 
low expected mass of the dark photon, the ideal tools 
to discover the dark photon and the dark Higgs boson 
are low energy and high-luminosity experiments such as 
Belle at KEKB, Belle II at SuperKEKB [ 23 ] , BaBar at 
PEP-II [25} [30]) and dedicated fixed target and beam 
dump experiments, several of which are planned or under 
construction [3llf37j . This article focuses on the Higgs- 
strahlung channel, e + e _ —i A'h'. Generally, the dark 
photon A' can decay into lepton pairs, hadrons, or in¬ 
visible particles while the dark Higgs boson h' can decay 
into either leptons pairs, or hadrons, where A'* 

is a virtual dark photon. The decay modes of the A' and 
h' depend on their masses and decay lengths [24, 381 . 
There are three main cases: (a) my < my: h' is 
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long-lived and decays to lepton pairs or hadrons, (b) 
to a' < mh' < 2mA''■ h' —► A'A'*, where A'* decays into 
leptons, and (c) > 2m a>'- h' —► A'A'. This article 

is concerned with case (c); in particular, we investigate 
ten exclusive final-states of type 3 2(1 + 1~)(tt + tt~), 

2(TT + n~)(l + l~), and 3(7r + 7r _ ), where l + l~ is an electron 
or muon pair but not a tau pair, and three inclusive final- 
states of type 2 (l + l~)X, where X is a dark photon can¬ 
didate detected via missing mass. 

The Higgs-strahlung channel involves the effective cou¬ 
pling of the dark photon to SM particles, a', induced via 
kinematic mixing with the SM photon, and the coupling 
of the dark-photon to the dark Higgs boson, an- KLOE 
and BaBar have reported searches for the dark photon 
and the dark Higgs boson [221 !2S : KLOE focused on 
mh' < mA' and BaBar on TOft,/ > 2mA' (assuming prompt 
decays of the A! and h r ), but no signal was found in ei¬ 
ther case. BaBar set limits on the product aj)X£ 2 (where 
e 2 = a'/a em and a e m is the SM electromagnetic coupling 
constant) for dark photon and dark Higgs boson mass 
ranges of 0.25 - 3.0 GeV/c 2 and 0.8 10.0 GeV/c 2 , re¬ 

spectively. Beam dump experiments [30H37) have placed 
90% confidence level upper limits on e for the processes 
e~p —> A'X' and pp —> A'X' (where X' is not identi¬ 
fied) of e < 10 -4 for a dark photon mass range of 1 - 
300 MeV/c 2 . Recently, BaBar [j3U] set an upper limit 
of e < 3 x 10~ 3 for a dark photon mass range of 0.3 - 
10 GeV/c 2 for the radiative decay process e + e _ —> yA'. 
The advantage of the Higgs-strahlung channel compared 
to the radiative decay is that the Quantum Electrody¬ 
namic (QED) background is expected to be much smaller. 
If, in addition, the coupling between the dark photon and 
the dark Higgs boson is of order unity, then the Higgs- 
strahlung channel is the most sensitive probe for the dark 
photon. 

Here, we report individual upper limits on the branch¬ 
ing fraction times the Born cross section, B x OBorin for 
the thirteen aforementioned Higgs-strahlung final states 
as well as combined upper limits on (JBorn and on the 
product an x e 2 for these final-states, in the mass ranges 
0.1 GeV/c 2 < toa' < 3.5 GeV/c 2 and 0.2 GeV/c 2 < 
mh' < 10.5 GeV/c 2 , assuming prompt decays of the dark 
particles. We use data collected with the Belle detec¬ 
tor [35] at the KERB e + e _ collider [35] . amounting to 
977 fhR 1 at center-of-mass energies corresponding to the 
T(1S) to T(5S) resonances and in the nearby continuum. 

We optimize the selection criteria and determine the 
e + e _ —> A'h' signal detection efficiency using a Monte 
Carlo (MC) simulation where the interaction kinematics 
and detector response are simulated with the packages 
MadGraph [50] and GEANT3 [51], respectively. There 
is no suitable background simulation available, so back¬ 
ground samples are taken from data sidebands. 

We choose loose particle identification criteria to en¬ 
hance the detection efficiency of final states with leptons. 
To ensure that only prompt decays are selected, i.e., that 



FIG. 1 : Example mb — m 3 ., distribution for the A'h! —» 

A cand A cand 

A'A'A' — » 67t channel, for m\, = 2.0 ± 0.1 GeV/c 2 , 

where m\, and m 3 A , are the dark photon candidates 

cand cand 

with the highest and lowest mass, respectively. The “same- 
sign” distributions (blue), where at least one A' candidate 
is reconstructed from 7r + 7r + or 7r _ 7r _ , are normalized to the 
“opposite-sign” 3(7r + 7r _ ) distributions (red) in the sidebands, 
and are used to predict the background in the signal region. 


the decay of each A' candidate occurs near the e + e _ in¬ 
teraction point (IP), we require that the vertex fit of all 
tracks detected in the event be consistent with an origin 
at the IP, and that each track have impact parameters 
\dz\ < 1.5 cm and dr < 0.2 cm, where dz is measured 
along the positron beam (collinear with the z axis) and 
dr is measured in the transverse r-(j> plane. We also re¬ 
quire that the second-order Fox-Wolfram moment [52] 
satisfy A 2 < 0.9, and that the electron helicity angle, a e , 
in the A' rest frame satisfy cos(a e ) < 0.9, as in Ref. |29j . 

For exclusive channels, we select final states with ex¬ 
actly three pairs of oppositely charged particles. For 
inclusive channels, we select final states of the type 
2{l + l~)X, where X is constrained by the missing mass 
of the event and contains zero, one, or two reconstructed 
tracks that are not identified as leptons or pions. We 
require that both TO;+;- and mx be greater than 1.1 
GeV/c 2 . Above this mass, the branching fraction of A! 
to hadronic final-states other than charged pion pairs is 
dominant [23] , We refer to events selected according to 
these criteria as “opposite-sign” to distinguish them from 
the “same-sign” events used for background estimation. 

For exclusive final states, we select candidate events 
with final-state masses between 98% and 105% of the 
initial-state mass. For inclusive channels, where this 
condition cannot be applied, we perform a missing-mass 
analysis: X is treated as an unobserved particle whose 
missing four-momentum is given by 

P >x = Pe+e~ ~ ^A' and ->Z + i- ~ + ’ W 

where P e + e - and P ^’, 2 _^ +J _ are the four-momenta of 
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FIG. 2: (a) : Signal candidates observed versus dark pho¬ 
ton candidate mass, m A i , and dark Higgs boson candidate 
mass, m A i A i , for the 13 final states. There are three 
entries per event, (b) and (c): Projection of signal candidates 
onto m A i A i and m A i (red points) with the predicted 
background (blue squares) from the scaled same-sign distri¬ 
butions for comparison. The dark photon candidate mass 
distribution has been scaled by 1/3. (d): Normalized resid¬ 
uals between the signal candidate distribution and predicted 
background, versus dark photon candidate mass (red points) 
and dark Higgs boson candidate mass (black squares). The 
same-sign error bars contain statistical and systematic errors. 
For empty bins, the systematic error is one event. 


the initial-state and the two fully reconstructed dark pho¬ 
ton candidates, respectively. The mass mx of the miss¬ 
ing four-momentum Px is then compared to the recon¬ 
structed masses of dark photon candidates 1 and 2 using: 

A m = m x -(rn\, l+l -+m\, l+l .)/2. (2) 

cand cand 

We select inclusive final-states by requiring 

Am m i n < A tti < Am max , (3) 

where the optimized limits Am m i n and A m max each de¬ 
pend on the measured mean mass of dark photon candi¬ 
dates 1 and 2 and on the particular final state. 

For exclusive (inclusive) final-states, we then require 
the invariant masses of dark photon candidates, m A i , 
to be consistent with three (two) distinct A! —> l + l~ or 
7r + 7r _ decays. Signal candidates with three (two) consis¬ 
tent dark photon masses are kept by requiring 

m T ad <m ^a„d <m A( and ’ 

where the optimized limits m™" and each de- 

cand cand 

pend on the measured mean mass of the three (two) fully 


reconstructed dark photon candidates and on the simu¬ 
lated width of the invariant mass distribution of the dark 
photon for that mass. 

For each event, if there is more than one signal can¬ 
didate that fulfills the selection criteria for a given final 
state, we select the candidate with the smallest Am. For 
exclusive channels, we use: Am = Sf Ato 2 with 

A mi = m\, ~(m\, +m\, + m 3 A , )/3. (5) 

cand cand cand cand 

If an event satisfies the selection criteria for multiple fi¬ 
nal states, we allocate the event to a single final state 
to ensure that the datasets for each final state are sta¬ 
tistically independent. This is accomplished by selecting 
the lowest numbered final-state category from the follow¬ 
ing list: (1) exclusive with 6 leptons, (2) exclusive with 
four leptons, (3) exclusive with two leptons, (4) exclusive 
with six pions, and (5) inclusive final-states. For the sig¬ 
nal MC simulation, the fraction of events with multiple 
signal candidates ranges from 7% to 15% in the chan¬ 
nels where we need to apply this ordering. For data, the 
fraction is below 0.5% in all final states. 

We optimize the event selection, including particle 
identification, the final-state mass requirements, and the 
parameters Arn m j n , Am max , m™ n and m A ? x using 

cand cand 

the signal MC simulation only. Events reconstructed as 
described above are used for signal. Background distri¬ 
butions are derived from the same event sample, by us¬ 
ing events where at least one dark photon candidate is 
reconstructed from two tracks with charges of the same 
sign, enforcing all selection criteria except charge con¬ 
servation. We refer to these as “same-sign” events. We 
verify that the background estimation is consistent with 
data as shown in Fig. [l] We generate MC with specific 
dark photon and dark Higgs boson masses and interpo¬ 
late between samples where necessary. The detection ef¬ 
ficiencies are 20% and 30%, on average, for the 3(e + e“) 
and 3(/i + /r _ ) final-states, respectively. 

For setting limits, we also estimate the background us¬ 
ing “same-sign” events, but in this case they are from 
experimental data. We sort the dark photon candidates 
by mass in descending order, m\, > m 2 A ; > m 3 A , , 

cand cand cand 

and calculate the mass difference m A , — m A , . We 

cand cand 

divide the data into different bins of m\, , with each 

cand 

bin analyzed separately. We divide the m\, — m 3 ,, 

. cand cand 

distribution into two regions: signal and sideband. The 
signal region size is determined by equation Q . The side¬ 
band region starts at 1.5 times and ends at 5.0 times the 
signal-region upper limit. Figure [l] shows the distribution 
of the mass difference m\, — m? A , for the bin m\, 

^cand c.-md rand 

= 2.0 ± 0.1 GeV/c 2 for the six-pion final-state. We as¬ 
sume that, in the absence of signal, the same-sign and the 
opposite-sign distributions have the same shape (but dif¬ 
ferent yields) in both the signal region and the sideband. 
Therefore, for each m A , bin, the same-sign distribu- 
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FIG. 3: Left: 90% CL upper limit on the product B x (JBorn for each of the 13 final states considered versus dark photon mass 
for different hypotheses for the dark Higgs boson mass. Black, red, green, blue and yellow curves correspond to rrih' = 1, 3, 5, 
7, and 9 GeV/c 2 , respectively, for exclusive channels and mh' = 3, 4, 5, 7, and 9 GeV/c 2 , respectively, for inclusive channels. 
Right: 90% CL upper limit on the cross section of e + e _ — > A'h', h' —> A'A' versus dark photon and dark Higgs boson mass. 


TABLE I: Number of events observed after all selection cri¬ 


teria are applied. 


Final-state 

Events 

Final-state 

Events 

3(e _ e + ) 

1 

2 ( M + / 0 ( e + e ~) 

1 

3 (P + P~) 

2 

2(/t + AO( 7r+7r_ ) 

1 

3(tv + tv~) 

147 

2(7r + 7r _ )(e + e _ ) 

5 

2{e + e~){p + p~) 

7 

2(n + n~)(p + p~) 

6 

2(e + e _ )(7r + 7r _ ) 

2 

( e + e~){p + p,~)(n + n~) 

7 

2(e+e-)A' 

572 

(e + e~)(p + p~)X 

30 

2(p+p~)X 

20 




tion (blue squares) is scaled so that the number of events 
in the sideband agree with the number of opposite-sign 
events (red points) in the sideband. The expected back¬ 
ground in the signal region is then the scaled number of 
events of the same-sign distribution in that region. This 
procedure is illustrated by Fig. [T] The opposite-sign and 
scaled same-sign distributions are consistent in the sig¬ 
nal region and the sideband. In the presence of signal, we 
would expect an excess of opposite-sign events over the 
predicted background in the signal region, as can be seen 
for the signal MC distribution. Figure [2] summarizes the 
background estimation. Figure [2] (a) shows the distribu¬ 
tion of events measured as a function of the dark photon 
candidate mass, , and the dark Higgs boson can¬ 

didate mass, rriA' A' ■ Table IT] shows the number of 
events observed after all selection criteria are applied. 

Figures[2](b) and (c) show the projections on the mass 
axis of the dark Higgs boson and dark photon, respec¬ 


tively. The number of events observed in the signal 
region, N a b s , and the number of predicted background 
events, IVbkgj are in good agreement. Their differences 
are quantified by the normalized residuals, shown in 
Fig. 2 (d) and defined as (iV obs - N hkg )/^jal hs + cr 2 kg , 
where cr 0 b s and Obkg are the standard deviations of the 
distributions. In all cases, the number of events observed 
is consistent with the background estimate. For exclu¬ 
sive final states, the background is mostly due to pro¬ 
cesses with p and w resonance particles, such as SM 2y 
processes. The discontinuity at 1.1 GeV/c 2 in Fig. [ 2 ] (c) 
is an artifact of the selection criteria. 

The upper limits on B x CTeom and (JBorn are calculated 
for ranges of rr ia' and mw , based on the signal MC mass 
resolution, with a Bayesian inference method with the 
use of Markov Chain Monte Carlo [53;. The number of 
observed events can be expressed as: 


N ohs = (JBorn • (1 + <5) • |1 — n| 2 ■ C ■ B ■ £ + N\ 


bkg, 


( 6 ) 


where 1 + S is an initial-state radiative correction fac¬ 
tor, |1 — n| 2 is the vacuum polarization factor, C is the 
luminosity, e is the detection efficiency, and N \>kg is the 
number of predicted background events. We calculate, 
for the exclusive (inclusive) channels, 1 + 6 using the for¬ 
mulae in Ref. HU and assuming the theoretical cross 
section is proportional to 1/s J24), where s is the square 
of the initial-state mass, and also assuming a cut-off value 
corresponding to 98% (a value between 20% and 90%) of 
the initial-state mass. 1 + 5 varies from 0.804 (0.93) to 
0.807 (1.17) depending on s and for the inclusive chan- 









































































6 



%o - 5 


= 0.3 GeV/c 2 


10 6 '—Belle upper limit 

10' 7 f - BaBar upper limit/ 
Belle expected 
limit 





10 7 r 
10‘ 8 r 

10 9 



m K = 5 GeV/c 2 

J oTTtTT 

m A , [GeV/c 2 ] 







FIG. 4: 90% CL upper limit on the product a_o x e 2 versus dark photon mass (top row) and dark Higgs boson mass (bottom 
row) for Belle (solid red curve) and BaBar [22] (dashed black curve). BaBar limits should be divided by (1 + 8) before being 
compared with Belle limits. The blue dotted curve, which coincides more or less with the solid red curve, shows the expected 
Belle limit. 


nels also the effective cut-off value. We use l+<5 = 0.8055 
(1.0) and include the variation as a systematic error in 
the upper limit calculation. The value of 11 — n(s)| 2 is 
taken from Ref. [551 [56] and varies between 0.9248 and 
1.072 depending on s. For B x <7Born and O’Born; loga¬ 
rithmic priors are used, and for 1 + <5, |1 — II| 2 , C, B , e, 
and iVbkg Gaussian priors are used to take into account 
the systematic uncertainty. In Fig. [3] the left panel shows 
the 90% credibility level (CL) 1 upper limits on B x cr B orn 
versus the dark photon mass, for different hypotheses of 
the dark Higgs boson mass, for each of the 13 final states 
considered, while the right panel shows the combined up¬ 
per limit on <7Bom for e + e~ — > Ah 1 versus the dark pho¬ 
ton and dark Higgs boson mass. For the combined limit, 
compared to BaBar, we use two extra channels, 3(7r + 7r - ) 
and 2(e + e~)X, which contribute 91% of our background. 
The branching fractions were taken from Ref. [24] . 

The combined limit can also be expressed as a limit 
on the product ap x e 2 by using the equations described 
in Ref. [24] . Figure [4] shows the 90% CL upper limits 
on ajj x e 2 for Belle, expected and measured, and for 
BaBar, for five different mass hypotheses for the dark 
Higgs boson (top row) and dark photon (bottom row) 
masses. Note that the BaBar limits were based on the 
visible cross section, rather than the Born cross section. 
For the expected limit, we assume: N 0 ^, s = -Wkg- 


1 In common High Energy Physics usage, this credibility level 

has been reported as “confidence level”, which is a frequentist- 

statistics term. 


The inclusion of the 3(7r + 7r _ ) final state dramatically 
improves the limit around the p and uj resonances. The 
dominant sources of systematic uncertainties are: the 
integrated luminosity (1%), branching fractions (4%), 
track identification (6%), particle identification efficiency 
(5%), detection efficiency (15%), background estimation 
(10%) and initial-state radiation (15%). All systematic 
uncertainties added in quadrature amount to 25%. 

In summary, we search for the dark photon and the 
dark Higgs boson in the mass ranges 0.1 3.5 GeV/c 2 

and 0.2 - 10.5 GeV/c 2 , respectively. No significant signal 
is observed. We obtain individual and combined 90% CL 
upper limits on the product of branching fraction times 
the Born cross section, B x OBorin on the Born cross sec¬ 
tion, (JBorn; and on the product of the dark photon cou¬ 
pling to the dark Higgs boson and the kinetic mixing be¬ 
tween the Standard Model photon and the dark photon, 
an x e 2 . These limits improve upon and cover wider mass 
ranges than previous experiments and the limits in the 
final-states 3(7r + 7r“) and 2(e + e“)X, where X is a dark 
photon candidate detected via missing mass, are the first 
limits placed by any experiment. For an equal to 1/137, 
rrih’ < 8 GeV/c 2 , and toa' < 1 GeV/c 2 , we exclude val¬ 
ues of the mixing parameter, e, above ~ 8 x 10” 4 . In 
the mass ranges, and for modes, where previous mea¬ 
surements from BaBar exist, the limits reported here are 
almost a factor of two smaller. Since the backgrounds are 
very low to non-existent, the improvement scales nearly 
linearly with the integrated luminosity. This bodes well 
for future searches with Belle II. 
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